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We consider the correlation of the electromagnetic field to determine spatial coherence inside a 
sonoluminescing bubble. We explicitly calculate the first order correlation function for two limiting 
cases of the excitation field: a blackbody spectrum and a discrete multifrequency spectrum. The 
correlation length for blackbody fields at temperatures between 3000 K and 10000 K is found to 
be on the order of the optical wavelength, increasing with decreasing temperatures. We predict 
spectral lines in the emission spectrum of single bubble sonoluminescence in cooler bubbles with 
interior temperatures below 10000 K. 

PACS numbers: 78.60.Mq, 42.50Fx, 34.80. Dp, 03.65. Bz 



Sonoluminescence is the outburst of a short light pulse 
from a noble gas bubble acoustically levitated in a liquid 
by ultrasound @,@]- Some of its salient features are the 
short duration of outburst with a pulscwidth — typically 
on the order of picoseconds the large number of 

photons per outburst — typically on the order of a mil- 
lion, and the wide emission spectrum over and probably 
beyond the visible range. Recently, it has been argued 
that cooperative decay of excited states of metastable 
molecules could explain the light emission so peculiar in 
SL §j§. 

The theory of collective decay or lasing of the 
population-inverted medium of molecules naturally con- 
tains the short timescale of the sonoluminescence pulse 
and the large number of photons emitted per outburst. In 
single bubble sonoluminescence, the emitted spectrum, 
however, is found to be a continuum, often devoid of 
any structure resembling the spectral lines of the gas 
molecules of either the solute or the solvent. It is nat- 
ural then to assume, though naively, that the emission 
mechanism is not compatible with a lasing process which 
is expected to be monochromatic. This argument is not 
correct for the reason that cooperative emission does in- 
deed occur for a collection of molecules with a broad dis- 
tribution of wavelengths. The input driving field which 
establishes this polychromatic distribution can be cor- 
related over a length comparable to the relevant scale 
in sonoluminescence. This scale defines the length over 
which molecules are correlated through the common elec- 
tromagnetic field. It also determines the number of cor- 
related molecules N e f / and hence the number of emitted 
photons. A large N e ff is essential in the proposed mech- 
anism |j]j7| as it determines the coherent part of the de- 
cay, Icoherent °c ^eff ' ^ ne incoherent part of the emission 
goes as N, the total number of emitting entities. 

In this paper we address this specific question per- 
taining to the excitation of the gas molecules necessary 
for population inversion. Considering that the gas in 
the sonoluminescing bubbles can reach a temperature of 
10000 K (~ 1 eV), or even higher it is important 



to determine if there is correlation in the electric field 
at these temperatures. Towards that end, we calculate 
the first order correlation function for a generalized ex- 
citation field in two limits of high (fc^T S> Tim) and low 
(ksT <C Tilo) interior temperature where the excitation 
is dominated by (a) a blackbody field and (b) a discrete 
multifrequency field arising from atomic/molecular lines, 
respectively. For the blackbody field, we find that the 
first-order correlation length in the electric field is on the 
order of the wavelength in the visible range, which allows 
the collection of excited molecules to be correlated over 
this length scale. With increasing temperature this corre- 
lation length is found to decrease substantially. Discrete 
multiple frequency fields corresponding to the spectral 
lines are also found to have long correlation lengths. We 
find that the correlation in the excitation is dominated 
by a blackbody field in hottor bubbles whereas discrete 
multiple frequency fields dominate at temperatures below 
10000 K. 

We predict that the spectral lines will be observable in 
single bubble sonoluminescence (SBSL) at temperatures 
below 10000 K; these discrete lines are expected to be 
prominent in the blue part of the emission spectrum. The 
discrete lines were not observed in the SBSL spectrum in 
previous experiments because of extremely high temper- 
atures inside the bubbles. This explains the experimen- 
tal observation [^T) that the sodium D-lines are visible 
in multiple bubble sonoluminescence (MBSL), which are 
now known to be at a typical gas temperature of 5000 K 
p^[ ; the D-lines were absent in the SBSL spectrum be- 
cause the single bubbles were much hotter. Our analysis 
allows us to understand, for the first time, the connec- 



tion between multiple bubble sonoluminescence [|12 13 
and single bubble sonoluminescence, namely they occur 
due to the same mechanism though in different temper- 
ature regimes. 

Coherence in an atom-field system is best understood 
by studying the electromagnetic field modes O-M. If 
N excited atoms distributed over a range of frequencies 
decay coherently, the photon number in the correspond- 
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ing modes of the number state of the electromagnetic 
field coupled to the atoms will change accordingly. Since 
the atoms and the field form a closed system, it is then 
equivalent and even useful to consider coherence in the 
field modes. 

The degree of coherence is often expressed by the cross 
correlation (E(ri,ti)E(r 2 ,t 2 )) and higher order correla- 
tions (E(ri,ti)E(r 2 ,t 2 ) ) in the field variable E(r,t). 

For example, the first order correlation is relevant for 
double-slit type interference experiments. The intensity 
is governed by the superposition of the amplitude con- 
tributions: E Q (t 2 -ti) = Ci.E(ri,ii) + c 2 E(r 2 ,t 2 ). This 
results in an additional part 2Re{ciE(ri,ti)c 2 E(r 2 ,t 2 )) 
to the intensity \Eg(t 2 — ti)\ 2 coming from first order 
interference. This interference term is a complex func- 
tion in general and it describes the first order coherence 
(two-point correlation function in the field). Similarly, 
higher order coherences are defined by the higher order 
correlation functions. For example, four-point correla- 
tion function defines the second order coherence relevant 
for the Hanbury-Brown and Twiss experiment. Our cal- 
culation of higher order coherence which signifies purely 
quantum mechanical effects will be presented elsewhere. 

The first order coherence function is defined by 



C(ri,r 2 ;ti,t 2 ) = {E i (r 1 ,t 1 )E j (r 2 ,t 2 )), 



(1) 



where (...) denotes the ensemble averaging through a 
density matrix opertaor p such that (EiEj) = Tr[pEiEj], 
and p — (1/Z)exp(—(3H). The definition of p is phys- 
ically important, and may include the active medium 
transitions. A multifrcquency electric field can be ex- 
panded in terms of the monochromatic fields Ei(r,t) = 
^12^ Eiu( r ,t)- If the temperature inside the bubble is 



very high, then the dominant contribution to the corre- 
lation function will come from the blackbody field, and 
hence, p could be determined for a known temperature 
T. At intermediate temperatures, the active medium be- 
comes important and the spectral contribution from the 
transitions in the given system must be taken into ac- 
count as well. In both cases p factorizes into operators 
for each frequency to and p = Yl u p u . Each u> component 
of the field may be assumed to be independent of each 
other, in which case the number state representation of 
the field reduces to product of the number states of in- 
dividual frequencies, i.e. |n Wl n W2 . . .) = ILjJra^). Then 
the correlation function has the following expression: 

{E^rut^E^M)) = 
^^(n 0J \E i u{r 1 ,t 1 )E ju (r 2 ,t 2 )\n 0J ){n 0J \p 0J \nu). (2) 

to n^j 

The product of fields at two points may be given as 

£?i(ri,ti)f? J -(r2,t2) = 
53 ^(a u 4e- iu(tl " fe) ^(ri)4(>-2) + h.c), (3) 



which follows from the definition of E(r,t) — 
l ^2 UJ uj(a ul exp[—iujt}A IJ j(r) + h.c.) assuming Coulomb 
gauge in the interaction picture. The diagonal terms in 
and aj, such as a^a^ are ignored in the following. 
The quantization for each mode u requires that a u] a\ J + 
a tj a u — 2Tj[2 n " + !]• Ignoring the second term, which 
represents vacuum fluctuations, one obtains an expres- 
sion for the symmetrized correlation function: 



J 



Tr[p(EiEj + EjEi)/*] = [ ^ exp[ 



kT 



kT 



}Re(exp[-iw(ti - t^A^v^A*^)). (4) 



In the above equation, for equal times, the summation 
over the occupation number n w yields the usual Planck's 
law. We consider the effect of active medium via a dis- 
tribution function: 



(5) 



where the strengths of different lines are given by g UJp 
and uo c is a cutoff frequency. This function will depend 
on the composition of the gas (solute) inside the bubble 
and the surrounding liquid (solvent). The effect of tem- 
perature and pressure inside the bubble can be taken into 
account by replacing the (5-function by a Lorentzian with 
a full width at half maximum corresponding to inhomo- 
geneous broadening. f(w,p) may also be a continuous 
function in case of an ionic spectrum which are realized 
in very hot bubbles. The final spectrum of emitted light 



would contain both f(u),p) and the Planck spectrum. At 
high temperatures, the Planck distribution would domi- 
nate over the emission line distribution. For the emission 
lines to be observable in SBSL, the temperature inside 
the bubble must satisfy the condition 

. hu>„ , . 

At high temperatures, the right hand side of the inequal- 
ity becomes equal to fc^T. In order to see the lines in the 
emission spectrum of SBSL at a wavelength A, the inte- 
rior temperature has to be reduced below T* = hc/Xks- 
The visible range of 200 nm — > 700 nm in wavelength 
corresponds to a temperature range of 71000 K — > 20000 
K. 

Eq. 4 gives the complete expression of correlation of 
fields in any arbitrary shape. For Ti = rj and ti — tj, it 
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reduces to the Planck distribution. The vector potential 
A(r) satisfies the Hclmholtz equation: 



V 2 A w +w 2 e(r)A w = 0, 



(7) 



with appropriate boundary conditions. For a gas bub- 
ble with an average dielectric constant ei immersed in a 
liquid with a dielectric constant £2- The solution of the 
Helmholtz equation is normally constructed from the lin- 
ear combination of the solutions of the scalar field equa- 
tion, requiring that Ej_, eE r , B r and Bj_//x be contin- 
uous across the bubble surface, i.e. at r = a. With 
the use of the boundary conditions, the transverse elec- 
tric (TE) and transverse magnetic (TM) mode solutions 
(A r , Ag, A < j,) TE - TM can be constructed inside and outside 
the bubble 

/. Blackbody field: Analytic evaluation of the coher- 
ence function with the exact solutions of Eq. (7) is rather 
complicated. Here we give results for the limiting case 
of wavelengths short compared to the minimum bubble 
radius, r>A. However, in standard SBSL experiments 
the minimum bubble radius is usually on the order of the 
wavelength. Even though r ~ A, closed-form expressions 
are possible only in this case. This condition is strictly 
satisfied in the case of millimeter-sized laser-induced bub- 
bles f|. 

Let us choose a direction for the wave vector k as 
(sin 9 cos <fi, sin 9 sm<fi, cos 9). The coherence function in 
the longitudinal BQ] direction is given by 



C[ on9 (r, t) 



1 



\3u> P 6(u -uj p )huj 



sinfcr 

TiLoduj [ — cos kr\ cos wr. (8) 

kr 



where lu — kc. If the thermal contribution dominates, 
T huj/kB, then the spectral distribution is essentially 
Planckian. This high temperature result was first ob- 
tained by Bourret pl| who derived the second-order elec- 
tric correlation tensor of blackbody radiation using tech- 
niques similar to those employed in the theory of isotropic 
turbulence of an incompressible fluid. The high temper- 
ature limit of the above result reduces to 



C[ on9 (r,0) 



1 



kdk , sin kr 



— 1 kr 



cos kr) 



oc 



:(1 



d I r 

or a a 



(9) 



where, a = crp, and t@ = Ti/kBT is the thermal time. 
L{x) — 7r/2 coth7rizi — \/2x is the standard Langevin func- 
tion. The expression in Eq. 9 is plotted in Fig. 1 for three 
different temperatures of the blackbody field. The scale 
over which correlation reduces to a characteristic value, 
the coherence length R^, is seen to vary between 0.5 to 
0.2 fim for bubble temperatures ranging from 3000 K to 
10000 K. This mesoscopic scale of coherence length is 
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FIG.l Variation of the normalized coherence function at 
three different temperatures. 

comparable to the compressed bubble size. At tempera- 
tures above 20000 K, the correlation length is found to 
be less than 0.1 /im. 

II. Discrete field: At low enough temperatures gfko p ^> 
kpT, coherence in the field of a set of discrete emission 
lines can be studied in a similar manner |23| . The coher- 



ence function C-[ 
equal to 

C l r 9 (k p ,r,0) 



long 



in the long wavelength limit becomes 



M! + o ( (V) 4 )]. (10) 



The r dependence of C l ° ng is parabolic, and a coher- 
ence length i?^ can be defined as the scale over which 
C l ° ng is reduced to a characteristic value, half of its max- 
imum. For a given wavelength A p , 



\p/ir. 



(11) 



In the lowest order, coherence is seen to be strongly 
dependent on the wavelength (A p = 27r/fc p ), 



long 



A 



-3 



(12) 



Hence SL from the discrete spectral lines in the blue part 
of the spectrum should be stronger than that in the red 
part. This is true for the coherent radiation in the SL 
emission spectrum. Observation of such a trend is a 
strong indication that the emitted light is indeed partly 
coherent. 

In the blackbody regime, the correlation length is much 
less than 0.05 /xm at a typical temperature of 40000 K; 
thus the effective number of molecules participating in 
collective decay reduces greatly as N e ff = p(47ri?^/3), 
where p is the average density. The sonoluminescence 
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emission spectrum is not dominated by coherent emis- 
sion. I incoherent 

~ N oc Rf ninl and I co herent ~ -^e/ / cx - 
i?^. For a typical minimum bubble radius R m in of 0.5 \im 
and a total A" of 10 6 , the incoherent emission is compa- 
rable to the coherent emission. At higher temperatures, 
incoherent emission dominates. 

Another corroborating evidence towards the transition 
from high temperature to low temperature behavior can 
be obtained from the temporal shape of the emitted light 
pulse. The pulse shape should progressively change from 
either sech 2 or oscillatory behavior (signifying ringing) 
to exponential decay which signifies incoherent emis- 

sionfjjlCjg. 

Let us summarize the three important aspects of the 
emission spectrum following our analysis: (a) Bubbles 
with an interior temperature much above 70000 K will 
typically display a blackbody emission spectrum in the 
entire visible range. As the temperature is reduced, spec- 
tral structure will appear in the blue part of the emis- 
sion spectrum, progressing towards the red part. Ex- 
perimentally, it is found that the spectrum is black- 
body and almost identical at long wavelengths, while it 
is gas specific at short wavelengths consistent with 
our picture. Below 10000 K, the emission spectrum will 
show atomic/molecular lines in the entire visible range, 
(b) The lines in the spectrum would be less strong in 
the red part as the correlation length for a discrete 
field decreases with increasing wavelength very sharply: 
Qiong ^ ( c ) Our analysis provides a natural con- 

nection between SBSL and MBSL surrounding the ques- 
tion as to why spectral lines are conspicuously absent in 
conventional SBSL experiments, but distinctly present in 
MBSL [ fil]]l2| . This is primarily because the gas temper- 
ature in MBSL is on the order of 5000 K, whereas in 
SBSL it's much higher. 

In conclusion, we show that the electromagnetic field 
inside a sonolumincscing bubble has a finite correlation 
even at high temperatures. At high temperatures (T 3> 
40000 K) the emission spectrum is expected to be pre- 
dominantly blackbody. At low temperatures (T < 10000 
K) discrete emission lines are predicted in single bubble 
sonoluminescence. Our analysis explains the parametric 
differences between multibubble sonoluminescence and 
single bubble sonoluminescence, exemplified by their very 
different emission spectra. 
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